Excited decuplet baryons from QCD sum rules 
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A calculation of the mass spectrum in the baryon decuplet sector is presented using the method 
of QCD sum rules. Sum rules are derived for states of spin-parity 3/2+ and 1/2+ using both 
the conventional method and a parity-projection method. The predictive ability of the sum rules is 
explored by a Monte-Carlo based analysis procedure in which the three phenomenological parameters 
(mass, coupling, threshold) are treated as free parameters and fitted simultaneously. Taken together, 
the results give an improved determination of the mass spectrum in this sector from the perspective 
of non-perturbative QCD. 
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I. INTRODUCTION 



A goal of hadronic physics is to understand the baryon 
spectrum from QCD, the underlying theory of the strong 
interaction. Experimentally, the push is fueled by the 
physics program at JLab, and other accelerator facilities. 
To this end, the time-honored QCD sum rule method jlj 
adds a useful theoretical tool. It is a non-perturbative 
approach to QCD that reveals a direct connection be- 
tween hadronic observables and the QCD vacuum struc- 
ture via a few universal parameters called vacuum con- 
densates (vacuum expectation values of QCD local oper- 
ators), the most important of which are the quark con- 
densate, the mixed condensate and the gluon conden- 
sate. The method is analytical, physically transparent, 
and has minimal model dependence with well- understood 
limitations inherent in the operator product expansion 
(OPE). It provides a complementary view of the same 
non-perturbative physics to the numerical approach of 
lattice QCD. The method was applied to the decuplet 
sector not long after it was introduced 0, 0, 0, IE El 
However, only limited attention has been paid to this 
sector since then. In Ref. jfl, a systematic study was 
made that incudes some A states. Later, progress came 
in the analysis of QCD sum rules by Monte-Carlo sam- 
pling of errors [j| . It was applied to the study of s pin-3/2 
branch of baryon decuplet |9| , magnetic moments |lCt lll| 
and other spin-3/2 excited baryons [12|. 

From a theoretical standpoint, the main issue is how 
to isolate sum rules that couple to a particular spin- 
parity, since the interpolating fields used to construct 
the spin-3/2 states contain both spin-3/2 and spin-1/2 
components, and couple to both parities. Separation of 
the spin-3/2 and spin-1/2 components can be achieved 
by examining the Dirac structures upon which the sum 
rules are based, as done in Ref. 0. However, the parities 
are still mixed. In the pole-plus-continuum ansatz of the 
method, isolation of a particular parity relies on strong 
cancellations in the excited part of the spectrum. With- 
out such cancellations, the sum rules usually suffer con- 
taminations from higher states, a problem that renders 
them unstable. A solution to separate the parities exactly 



was proposed in Ref. |T^ which showed improved conver- 
gence in the octet baryon sector. In this work, we apply 
the method to the decuplet sector, together with the con- 
ventional method. The goal is to identify sum rules that 
couple strongly to the four branches of the decuplet with 
spin-parity 3/2+, 3/2—, 1/2+, and 1/2—. In addition 
to the Delta states, we also examine states that contain 
the strange quark, which have not been studied in detail. 
In our calculation, we consistently include operators up 
to dimension eight, first order strange quark mass cor- 
rections, flavor symmetry breaking of the strange quark 
condensates, anomalous dimension corrections, and pos- 
sible factorization violation of the four-quark condensate. 
Furthermore, we try to assess quantitatively the errors 
in the phenomenological parameters, using the Monte- 
Carlo based analysis procedure. This procedure incor- 
porates all uncertainties in the QCD input parameters 
simultaneously, and translates them into uncertainties in 
the phenomenological parameters, with careful regard to 
OPE convergence and ground state dominance. 



II. METHOD 

In the conventional method, the starting point is the 
time-ordered two-point correlation function in the QCD 
vacuum 

II(p) = i J d A x e lp - x (0 | T{ rj(x) fj(Q) } | 0). (1) 

where rj is the interpolating field with the quantum num- 
bers of the baryon under consideration. Assuming SU(2) 
isospin symmetry in the u and d quarks, we consider the 
most general current of spin 3/2 and isospin 3/2 for the 
A, 
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(x)Ca fl u b {x)]u c (x) 



(2) 



Here C is the charge conjugation operator and the su- 
perscript T means transpose. The antisymmetric e abc 
and sum over color ensures a color-singlet state. For the 
other members of the decuplet, we consider, omitting the 
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explicit ir-dependence 



v f = y/T/3e abc {2[u aT Ca^s b }u c + [u aT C a , (3) 



vf = VI73e Qbc {2[ S QT CaX].s c + [s aT Ca^s b ]u c }, (4) 



e abc \s aT C<7„s b ]s c 



(5) 



A baryon interpolating field couples to both the ground 
state and the excited states of a baryon, and to both par- 
ities. That leads to sum rules that contain mixed par- 
ity states. The parities can be separated by considering 
the 'forward-propagating' version of the correlation func- 
tion H3, 

n(p) =i I d 4 x e ip x 9(x ){Q \ T{ r,{x) fj(0) } | 0). (6) 



where the only difference between the above equation 
and the conventional two-point correlation function in 
Eq. (JTJ is the step function 0(xq). Under this condition, 
the phenomenological representation of the the imaginary 
part in the rest frame (p = 0) can be written as 



Imn(p ) = E 
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-S(p a - m+) 



2 7o 



(A") 



loA(p Q ) + B( Po ) 



■5(po-m n ) 



(7) 



where A 1 * 1 and are the couplings and masses of the 
states involved, separated by parity. The functions in the 
last step are defined by 



and 



(8) 



B M = \ E t( A « - m n) - {K?5{po m~)]. 

n 

(9) 

By considering the combinations A(po) + B(po) and 
A(po) — B(po), positive-parity and negative-parity states 
are separated exactly, respectively. This method is rem- 
iniscent of the one used in lattice QCD under Dirichlct 
boundary conditions in the time direction |l4j |. 

The construction of sum rules in Borel space proceeds 
by taking the integral with a weighting factor and trun- 
cating all excited states starting at a certain threshold. 
For positive-parity states, 

r w + 

/ [A( Po ) + B(p Q )} e-vH M * d Po = \\er m V M2 (10) 



and for positive-parity states, 

" [A( Po ) - B( Po )] e-P»/ M2 d Po = X 2 _e- m -/ M ' 2 . (11) 



To separate out the spin components, we need the full 
Dirac structure of the correlation function for spin-3/2 
interpolating fields 0, @ , 
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(12) 



where the ellipses denote excited state contributions, and 
p denotes p M 7^. The upper/lower sign corresponds to 
the explicit positive/negative parity when a state con- 
tributes. The entire structure is divided into two parts: 
one for spin-3/2 with coupling A^ 2 , and the other spin- 

1/2 with coupling a^ 2 . The structures are chosen to 
be effectively orthogonal to each other. Only six of the 
structures are independent. The two sum rules from g^ v p 
and g^ u couple only to spin-3/2 states because they do 
not appear in the spin- 1/2 part. The two sum rules from 
ItiPu + IvPv and 7 M 7„ + \g^ y couple only to spin- 1/2 
states because they do not appear in the spin-3/2 part. 
By the same token, the sum rule at 7^p„ + ^vPp, that 
couples to spin-1/2 could also be obtained from the com- 
bination 7 M 7„p + \giivP- The two sum rules from p^p^p 
and p^Pu couple to both spin-3/2 and spin-1/2 states 
because they appear in both parts. They will not be 
considered further. Note that the parities are still mixed 
in each of these sum rules. It is not known a priori which 
parity is the dominant state in a given sum rule. That 
information comes from the specific OPE structure on 
the left-hand-side. 



III. PARITY-PROJECTED QCD SUM RULES 
FOR SPIN-3/2 DECUPLET STATES 

The combination of the two tensor structures that cou- 
ple purely to spin-3/2 states, g^P and g M „, can be cast 
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into the form g^ u (joA + B) in the rest frame. According 
to Eq. 0, it suggests that parity projection can be done 
exactly in this case. The functions A(p ) and B(p ) are 
readily identified from the usual calculation of the OPE 
which involves contracting out the quark pairs in the cor- 
relation function and substituting the fully-interacting 
quark propagator from OPE. They are functions of the 
QCD vacuum condensates and other QCD parameters. 

Upon Borel integration in Eq. (|10|l and Eq. 1|11|) . the 
parity-projected sum rules can be written in the general 
form 

A(w+,M) + B(w+,M) = \ 2 + e- m + /M \ (13) 

A{w-,M) - B(w-,M) = X 2 _e- m -' M \ (14) 

where \± — (27r) 2 A+_ are the rescaled couplings. The 
rescaling is done so that no factors of tt appear explic- 
itly in the sum rules. The rescaled quantities are also 
'natural' in their numerical values. 

For our special case of spin-3/2 states, the A function 
A 3 / 2 is given by 

A 3/2 = c 1 [l](2(l-e- tu2 / M2 )M 6 -2w 2 e- w2 ' M2 M i 

- W 4 e- W ^ M2 M 2 )L^ 27 

+ (ca [b] + c 3 [m s a]) (1 - e— 2 /M 3 )M 2 L 4/27 

+ c 4 [m s m 2 a] L^ 27 + c 5 [n v a 2 } L 28 / 27 . (15) 

In this expression, the vacuum condensates are explicitly 
isolated in square brackets whose definitions and values 
will be given below. The Wilson coefficients (which are 
pure numbers) are given for various members of the de- 
cuplet by 
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Similarly, the B 3 / 2 function is given by 

B 3/2 = ci[m s ][I(w) -we~ w2 / M2 )L- 8 / 27 M 4 

+ c 2 {m s }w 3 e- w2 / M2 L~ 8 / 27 M 2 

+ c 3 [a](I(w) - we- w2 ' M2 )L w ' 27 M 2 

+ c 4 [m s b]I{w)L- 8 ' 27 

+ c b [m 2 a]I{w)L 2 ' 27 ', (17) 



where the coefficients are given by 
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The integral I(w) — J™ e x l M dx is evaluated numer- 
ically in the analysis. Note that the sum rule for Q,~ 
reduces to that for the A if m s = and f s = 1, which 
serves as a check of the calculation. 

Now we explain the meaning of the parameters in the 
sum rules. The rescaled quark condensate is taken as 
the positive quantity a = — (2tt) 2 (uu) = 0.52 ± 0.05 
GeV 3 , corresponding to a central value of (uu) = — (236) 3 
MeV 3 . For the gluon condensate, we use rescaled b = 
{g 2 c G 2 ) = 1.2 ± 0.6 GeV 4 . The mixed condensate pa- 
rameter is placed at Wq = (qga ■ Gq)/(qq) = 0.72 ± 0.08 
GeV 2 . For the four-quark condensate {uuuu) = k v (uu) 2 , 
we use k v = 2 ± 1 to include possible violation of the fac- 
torization approximation. We retain only terms linear 
in the strange quark mass and set m u = = 0. The 
strange quark mass is taken as m s = 0.15±0.02 GeV. We 
use the ratio f s = (ss)/(uu) = (sg c a ■ Gs) / {ug c a ■ Gu) = 
0.83 ± 0.05 to accounts for the flavor symmetry break- 
ing of the strange quark condensates. The anomalous 
dimension corrections of the various operators are taken 
into account via the factors L 7 = [a s (/i 2 )/a! s (M 2 )] 7 = 
[ln(M 2 /A 2 3CD )/ln(^ 2 /A 2 3CI) )] 7 , where 7 is the appro- 
priate anomalous dimension, /1 = 500 MeV is the renor- 
malization scale, and Aqcd = 0.15 ± 0.04 GeV is the 
QCD scale parameter. We find variations of Aqcd have 
little effects on the results. The uncertainties are assigned 
fairly conservatively, ranging from 10% to 100%. They 
will be mapped into those for the fit parameters in the 
analysis, giving a realistic estimate of the errors on the 
parameters. 

It is worth pointing out that the A function is chiral- 
even, meaning that it only involves condensates of even 
energy dimensions: for example, & of 4, m s a of 4, n v a 2 
of 6 and so on. On the other hand, the B function is 
chiral-odd, with the leading contribution from the quark 
condensate a, the order parameter of spontaneous chiral 
symmetry breaking of QCD. The difference of the two 
parity-projected sum rules in Eq. (|T3^1 and Eq. ifHI) is 
the sign in front of the chiral-odd term B. So from the 
point of view of QCD sum rules, the origin of splittings 
between positive-parity and negative-parity states lies in 
the chiral-odd vacuum condensates. If the quark con- 
densate and the strange quark mass vanish, then chiral 
symmetry is restored in the vacuum, and there would be 
exact parity doubling in the baryon spectrum. This is 
a valuable insight from the parity-projected QCD sum 
rules. 

The analysis of a QCD sum rule boils down to the 
following mathematical problem. Given an equation of 
the general structure 

LHS(Af, OPE) = RHS(M, w, m, A 2 ), (19) 

and a set of input QCD parameters denoted by OPE, 
find the best output parameters (the baryon mass of in- 
terest 77i, the coupling strength A 2 of the interpolating 
field, and the continuum threshold w) by matching the 
two sides over a region in Borel mass M. The LHS has 
errors arising from our imprecise knowledge of the QCD 
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parameters, as discussed above. From statistical point of 
view, a x 2 minimization of the type 

9 x - ^ |LHSj — RHSjl 2 .„„. 
X 2 = E J —2 -« ( 2 °) 

offers the least-biased way of finding the unknown pa- 
rameters. In the Monte-Carlo based procedure used here, 
first the Borel window in M is divided equally into a grid 
(we use 51 points). At each point, the uncertainty distri- 
bution in the OPE is constructed by randomly-selected, 
Gaussianly-distributed sets generated from the central 
values and the uncertainties in the QCD input param- 
eters. Then the x 2 minimization is applied to the sum 
rule by fitting the phenomenological parameters. This is 
done for each QCD parameter set and at each point in 
M, resulting in distributions for phenomenological fit pa- 
rameters, from which their errors are derived. Usually, 
100 such configurations are sufficient for getting stable 
results. We generally select 500 sets which help resolve 
more subtle correlations among the QCD parameters and 
the phenomenological fit parameters. So the same sum 
rule is fitted not just once, but thousands of times in our 
analysis. The advantage afforded by the Monte-Carlo is 
that the entire phase-space of the input QCD parame- 
ters is explored. The errors obtained this way represent 
the most realistic and conservative estimates of the pre- 
dictive power of a QCD sum rule. This is in contrast 
to traditional approaches where only a small part of the 
phase-space is explored at a time. We select the Borel 
window by trial and error. Our criteria are that the OPE 
is reasonably convergent by looking at the terms of var- 
ious dimension, and that the fit results should not be 
sensitive to small changes in the Borel window. Under 
these constraints, we seek to make the window as wide as 
possible. Since the selection of Borel window M and the 
three phenomenological parameters are inter-dependent 
in the sum rule, the entire fitting process is iterated until 
the best solution is found. In general, we seek solutions 
that are from a three-parameter search in which all three 
parameters are treated as free simultaneously. We regard 
such solutions as the best predictions of the QCD sum 
rule approach. In the absence of such solutions, a two- 
parameter search is performed in which one of the three 
parameters is fixed. Usually, it is the continuum thresh- 
old which is fixed to values that are larger than the mass 
in question, or to values suggested by the observed spec- 
trum. Such a two-parameter approach is the one usually 
adopted in most analyses in the past. 

Fig.nshows the matching in the sum rule Eq. (|13|) for 
3/2+ states. The two sides match very well over a rela- 
tively wide region (1 GeV). The individual contributions 
of term A and term B are also shown. The two terms 
are comparable in size, the sum of which gives rise to the 
LHS. The error band on the LHS is generated by Monte- 
Carlo reflecting all the uncertainties assigned to the QCD 
parameters. The entire uncertainty phase space is then 
mapped to the output parameters in the fitting process. 
A three-parameter search was tried first. Unfortunately, 



TABLE I: Results for the 3/2+ states from the parity- 
projected sum rule Eq.( I13P in a two-parameter search. The 
errors are derived from 500 Monte-Carlo sets based on the un- 
certainties assigned to the QCD input parameters. Two sets 
of solutions are provided for each case: the first one using the 
default errors in the OPE input parameters, while the second 
using uniform 10% errors. The experimental values are taken 
from the PDG 





Region w -^3/2+ Mass Exp. 
(GeV) (GeV) (GeV 6 ) (GeV) (GeV) 


A(|+) 


1.4 to 2.4 1.60 2.46+0.42 1.21+0.06 1.23 
2.32+0.21 1.21+0.02 


E*(§+) 


1.4 to 2.4 1.84 3.52+0.43 1.38+0.06 1.385 
3.40+0.29 1.37+0.02 


-*(§+) 


1.4 to 2.4 1.95 4.19+0.46 1.47+0.05 1.53 
4.10+0.37 1.47+0.02 


«-(§+) 


1.4 to 2.4 2.25 6.94+0.63 1.68+0.05 1.67 
6.89+0.56 1.68+0.02 



a solution could not be found. What happens is that the 
search algorithm keeps returning solutions with thresh- 
old smaller than the mass, a clearly unphysical situa- 
tion. This is an indication that the OPE does not have 
enough information to resolve all three parameters simul- 
taneously. So we switched to a two-parameter search by 
fixing the continuum threshold to values suggested by the 
Particle Data Group ^5|. The extracted parameters are 
given in Table |U We offer two sets of results on the fit 
parameters: one using the default errors assigned to the 
OPE parameters ranging from 10% to 100%, the other 
with 10% uniform errors. It sort of gives the worst-case, 
best-case scenarios, as far as errors on the input parame- 
ters are concerned. In the worst-case solution, the errors 
on the masses are on the order of 5%, and 15% on the 
couplings. In the best-case solution, the errors are re- 
duced to about 2% on the masses, and less than 10% 
on the couplings. Note that the central value for the 
couplings are shifted by a small amount when the errors 
are reduced, while it is stable for the masses. The com- 
puted masses compare favorably with experiment. What 
is interesting is the fact that the mass pattern emerges 
under the same Borel window across the particles, and 
that the continuum thresholds are consistent with the ex- 
cited states suggested by PDG. In the case of S*(3/+), 
the computed mass of 1.47 GeV is slightly below the ex- 
perimental value of 1.53 GeV. In the PDG, there are two 
3-star excited states at 1690 MeV and 1950 MeV with 
unknown spin-parity. Our sum rule favors a threshold of 
1950 MeV to 1690 MeV which leads to a mass that is 
even smaller than 1.47 GeV. It hints that the spin-parity 
of 1690 MeV is likely not 3/2+. 

Fig.|21shows the matching in the sum rule Eq. (|14|) for 
3/2— states. Here it is the difference between term A 
and term B that gives rise to the ground-state pole. The 
quality of the matching is still fairly good, but we found 
that the stability of the sum rules is not as good as the 
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FIG. 1: Matching of the sum rules in Eq. IUM|I for positive-parity spin-3/2 states as a function of the Borel mass. The solid line 
along with the error band generated via Monte-Carlo is the LHS (A + £>)• The dash-dot line, which is barely distinguishable 
from the solid line, is the RHS (ground-state pole). Also plotted are the individual contributions of term A (dotted line) and 
term B (dashed line). 



corresponding ones for 3/2+ states in Eq. i|13|) . A three- 
parameter search does not work, so a two-parameter 
search is performed. To get a better understanding of 
the sum rules, we adopt the following fitting strategy. 
Instead of fixing the continuum threshold to a value 
that leads to a known state, we perform two-parameter 
searches by fixing the mass to values suggested by the 
PDG, and search for the continuum threshold and cou- 
pling. In this way, the searches will provide informa- 
tion on whether a specific sum rule can accommodate a 
known state with reasonable continuum thresholds and 
couplings. Such a study is useful since to our knowledge 
there is no information about these 3/2— states from the 
standpoint of QCD sum rules. The result of the study 
is given in Table ITT1 We fix the Borel window to be the 
same wide window (1.4 to 2.4 GeV) as for the 3/2+ states 
for all cases. This choice is quite reasonable judging by 
the matching plots, which are all of comparable quality 
as Fig. |3 Again, the same worst-case, best-case solu- 
tions are provided for each case in order to give some 



idea about the stability of the fits. 

In the case of A(3/2— ), the PDG lists a 4-star state at 
1.70 GeV. Using this value as input, a continuum thresh- 
old of 2.74 GeV and a coupling of 1.56 are obtained, with 
relative errors about 20% and 55%, respectively. They 
decrease to about 5% and 40% in the best-case (10% 
uniform errors on the input parameters). Note that the 
central values vary as errors are reduced: the threshold 
shifts up to 2.9 GeV and the coupling shifts down to 
1.43. It is a sign of highly non-linear mapping of the er- 
rors from input to output. These errors are much bigger 
than those for 3/2+ states, especially on the couplings. 
The next A(3/2— ) state in the PDG is a one-star state 
at 1.94 GeV, compared to the continuum threshold of 2.7 
to 2.9 GeV from this sum rule. 

In the case of E*(3/2— ), the PDG lists a 2-star state 
at 1.58 GeV with unknown spin-parity, a 4-star state 
at 1.67 GeV, and a three-tar state at 1.94 GeV. The 
corresponding continuum thresholds required to predict 
these states are 2.87, 2.90, 3.00 GeV. 
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TABLE II: Results for the 3/2— states from the parity- 
projected sum rule Eq. ( I14p in a two-parameter search. The 
errors are derived from 500 Monte-Carlo sets based on the un- 
certainties assigned to the QCD input parameters. Two sets 
of solutions are provided for each case: the first one using the 
default errors in the OPE input parameters, while the second 
using uniform 10% errors. The experimental values are taken 
from the PDG El- 





Region 




^3/2— 


Mass 


Exp. 




y\j^ v ) 


v v ) 


V ) 


(GeV) (GeV) 




1 4 tn 9 A 


9 74+n A7 

9 Qn+n 1 n 


1 ^fi+n 84 
1 4^+n ^fi 


1 TC\ 
1. / U 


1. ( U 


y*(3 \ 


1 A 9 A 


9 fis+n m 
2.87±0.13 


1 1 ^+n fi4 
1.02±0.39 


l.Do 








9 74+n A(\ 


1 ^1 +n 79 


l.o ( 


l.D( 






9 Qn+n 1 1 


1 1 Q+n 4fi 










9 Q4+n ^n 


9 14+1 on 


1 Q/1 








3 nn+n ofi 


1 Q4+n 7fi 






•B*(l—\ 
- \2 > 


1.4 to 2.4 


79+n 4*3 
2.91±0.12 


1.02±0.42 


1.69 


1.69(?) 






2.86±0.36 


1.40±0.81 


1.82 


1.82 






2.95±0.10 


1.28±0.53 










2.94±0.30 


1.80±1.05 


1.95 


1.95(?) 






3.00±0.10 


1.63±0.68 






fi-(§-) 


1.4 to 2.4 


3.11±0.13 


2.31±1.11 


2.25 


2.25(?) 






3.18±0.24 


3.47±2.27 


2.38 


2.38(?) 






3.18±0.13 


3.13±1.50 










3.24±0.24 


4.32±2.77 


2.47 


2.47(?) 






3.24±0.14 


3.91±1.82 







In the case of S*(3/2— ), out of the three 3-star states 
listed by the PDG, 1.69, 1.82, 1.95 GeV, only the 1.82 
GeV state is assigned the 3/2— spin-parity. The corre- 
sponding continuum thresholds required to predict these 
states are 2.91, 2.95, 3.00 GeV. Interestingly, these values 
are close to the ones for the £*(3/2— ) states. 

In the case of f2~(3/2— ), there are three candidates 
from the PDG, a three-star state at 2.25 GeV, a two- 
star state at 2.38 GeV, and a two-star state at 2.47 GeV, 
all with unknown spin-parity. Using the sum rules, they 
could be assigned 3/2— if the the corresponding contin- 
uum thresholds are 3.11, 3.18, 3.24 GeV. For the 2.25 
GeV state, we could find a solution only when the errors 
are small. 



IV. CONVENTIONAL QCD SUM RULES FOR 
SPIN-1/2 DECUPLET STATES 

The two sum rules that couple purely to spin- 1/2 
states are from the Dirac structures 7 M p„ + IvVix and 
l^liy + \g^v Since they cannot be cast into the (joA+B) 
form in the rest frame, the parity projection technique 



cannot be applied to this case. We have to rely on 
the conventional sum rule method in which the parities 
are mixed. The spin-1/2 decuplets are in fact excited 
states. They are rarely studied in the QCD sum rule 
method. Here we have an opportunity to isolate them as 
the ground-state poles. For this reason, even a conven- 
tional analysis is beneficial. 

The sum rule from the ^^Pu + IvPn structure is 



Cl [1] L 4 / 27 E 2 M 6 
+c 3 [m s a 



c 2 [b] L 4 / 27 E Q M 2 



L A ' 27 E M 2 + c 4 [m s mla] L" 10 / 27 



+c 5 [K v a 



"1/2- 



l 28 / 27 + C6 [ m y 

where the coefficients are given by 



r 14/27 1 

M 2 

-m 2 1/2+ /M 2 



(21) 
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— — 
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216 
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R48 



7/ a (l+2/ 3 ) 



111 

21fi 



This sum rule is chiral-even. Since the parity of the low- 
est state is not known, we retain one term for each parity 
in the RHS, and let the OPE reveal which one is lower. 
The sum rule for A agrees with that given in Ref. 0, 
except that the coefficient cq has the opposite sign. The 
sum rules for £*, S *, and f2~ are new, to the best of 
our knowledge. Note that the sum rule for fl~ reduces 
to that for A if the strange quark is turned off (m s = 
and f s = 1), as expected. Since they are derived sepa- 
rately, this provides a non-trial check of the calculation. 
Another check is provided by the fact that C5 and Cq for 
£*, S*, and coincide with each other in the limit of 
f s = 1, as expected. The excited state contributions of 
RHS are modeled using terms on the OPE side surviving 
AI 2 — » 00 under the assumption of duality, and are rep- 
resented by the factors E n {x) = 1 — e~ x ^ x n /n\ with 
x = w 2 /M 2 and w an effective continuum threshold. 

The sum rule from the 7^7^ + \g^i/ structure, which is 
chiral-odd, is given by 

ci [a] L 16/27 E x M 4 + c 2 [ 
[m s b] L- 8 ' 27 E M 2 

^/2+/M 2 



-C3 

-C5 



c 4 



] I 2 ' 27 E M 2 
[ab] L 16 / 27 



-& 2 



l/ 2+ '«l/2+t 



+a 2 



'1/2 



_/M 2 



l/2- m l/2- 

where the coefficients are given by 



(23) 
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(24) 
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FIG. 2: Matching of the sum rules in Eq. il l 1ft for negative-parity spin-3/2 states. The solid line along with the error band 
generated via Monte-Carlo is the LHS (.4 — 23). The dash-dot line, which is barely distinguishable from the solid line, is the 
RHS (ground-state pole). Also plotted are the individual contributions of term A (dotted line) and term B (dashed line). 



For A, the coefficient C4 is different in both value and 
sign from that in Ref. 0. The sum rules for £*, S*, 
and ft~ are new, as far as we know. Similarly, the sum 
rule for f2 _ reduces to that for A if the strange quark is 
turned off. Another check is provided by the fact that 
c\, C2, C4 and C5 for £*, S*, and Or coincide with each 
other in the limit of f s = 1. 

Now we turn to the analysis of the two sum rules. 
First, we note that states of opposite parities on the phe- 
nomenological side (RHS) are adding up in the chiral- 
even sum rule, whereas canceling in the chiral-odd sum 
rule. This is a standard feature of baryon sum rules 
that leads to the general conclusion that chiral-odd sum 
rules perform better than chiral-even sum rules in baryon 
channels. Indeed, we found that the chiral-even sum rule 
in Eq. I|21|) is very poor. The leading term is a perturba- 
tive contribution, and the sum rule is almost completely 
saturated by the continuum. No results could be ex- 
tracted from this sum rule. 

On the other hand, the chiral-odd sum rule in Eq. i|23l) 
has good convergence. The leading term contains the 



non-perturbative quark condensate. The sign of this 
term (in ci) indicates that the sum rule is saturated by 
the positive-parity state: they have the same negative 
signs, see Eq. <\2'.il and Eq. J23J. This is an example of 
how the parity of the ground-state pole is determined in 
a mixed-parity sum rule. The dominance of the quark- 
condensate term is further confirmed in our numerical 
Monte-Carlo analysis. We are able to perform three- 
parameter searches. It means that the extracted mass, 
coupling, continuum threshold can be regarded as the 
true predictions of the sum rule. The results are given in 
Table Hill and the matching of the two sides is shown in 
Fig- El lu the observed spectrum, the lowest A(l/2+) is 
a one-star state at 1.75 GeV, followed by a 4-star state at 
1.91 GeV. Our prediction favors the latter as the lowest 
state in this channel. It casts doubts on the existence of 
the state at 1.75 GeV. For E*(l/2+), the PDG lists a 
3-star state at 1.66 GeV, followed by a one-star state at 
1.77 GeV, then by a 2-star state at 1.88 GeV. Our re- 
sult seems to favor the state at 1.88 GeV than the state 
at 1.66 GeV. The spin-parity situation in the H*(l/2+) 
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TABLE III: Results for the 1/2+ branch of the decuplet states 
from the chiral-odd sum rule Eq. J23 

in a three-parameter 
search. The errors are derived from 500 Monte-Carlo sets 
based on the uncertainties assigned to the QCD input pa- 
rameters. Two sets of solutions are provided for each case: 
the first one using the default errors in the OPE input pa- 
rameters, while the second using uniform 10% errors. The 
experimental values are taken from the PDG |l5j|. 

Region w &1/2+ Mass Exp. 
(GeV) (GeV) (GeV 6 ) (GeV) (GeV) 

A(l/2+) 1.1 to 2.0 2.41±1.08 0.18±0.15 2.10±0.54 1.91 
2.52+0.80 0.18+0.13 2.17+0.33 

E*(l/2+) 1.1 to 2.0 2.15+1.06 0.13+0.13 1.97+0.62 1.88 
2.29+0.83 0.13+0.11 2.09+0.37 

S*(l/2+) 1.1 to 2.0 2.42+1.11 0.16+0.13 2.11+0.54 1.95(?) 
2.50+1.50 0.16+0.12 2.19+0.35 

fr(l/2+) 1.1 to 2.0 3.02+1.18 0.23+0.14 2.37+0.41 2.38(?) 

3.11+1.13 0.23+0.13 2.39+0.32 



channel is not clear in the PDG. It lists two 3-star states 
with unknown spin-parity, at 1690 MeV and 1950 MeV. 
Our results is in favor of the 1950 MeV state. In the f2~ 
channel, there are 3 states sitting close to each other in 
the observed spectrum, at 2250, 2380, and 2470 MeV, 
whose spin and parity are not clear. Our prediction of 
2.37 GeV is in the middle of this range, but the accuracy 
is not enough to clearly identify with one of them. 



V. CONCLUSION 

We have presented a study of the decuplet family using 
the method of QCD sum rules. New sum rules are derived 
for the 3/2+ and 3/2— branches using a parity-projection 
technique. They are more stable and give a better de- 
termination of the mass spectrum than the conventional 
sum rules. The results for the 3/2— branch are new as a 
consequence of the parity separation. The spin-1/2 sector 
is investigated using the conventional sum rules method. 
New sum rules are derived for members that contain the 
strange quark (E*, 3*, and in addition to a care- 

ful re-examination of A channel. The chiral-even sum 
rules of Eq. 121|) are dominated by the continuum. No 
useful information could be extracted from them. The 
chiral-odd sum rules of Eq. &2'A\\ have good convergence 
and allow the only three-parameter searches in this study. 
The predicted results provide useful information on the 
1/2+ states from a QCD-based standpoint. We could 
not find sum rules that are saturated by the 1/2— states. 
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FIG. 3: Matching of the sum rules in Eq. i'2'-'>t for the spin-1/2 decuplet states. The solid or black line along with the error 
band generated via Monte-Carlo is the LHS (OPE minus continuum). The dash-dot or red line, which is barely distinguishable 
from the solid line, is the RHS (ground-state pole). Also plotted are the individual contributions of the OPE (dashed or green 
line) and the continuum (dotted or blue line). 



